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In this paper, the dynamical equations and junction conditions at the interface between adjacent 
layers of different elastic properties for an elastic deformable astronomical body in the first post- 
Newtonian approximation of Einstein theory of gravity are discussed in both rotating Cartesian 
coordinates and rotating spherical coordinates. The unperturbed rotating body (the ground state) 
is described as uniformly rotating, stationary and axisymmetric configuration in an asymptotically 
flat space-time manifold. Deviations from the equilibrium configuration are described by means of a 
displacement field. In terms of the formalism of relativistic celestial mechanics developed by Damour, 
Soffel & Xu, and the framework established by Carter & Quintana the post Newtonian equations 
of the displacement field and the symmetric trace-free shear tensor are obtained. Corresponding 
post-Newtonian junction conditions at interfaces also the outer surface boundary conditions are 
presented. The PN junction condition is an extension of Wahr's one which is a Newtonian junction 
conditions without rotating. 
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I. INTRODUCTION 

The theory of elastic deformable bodies is of great importance for quantitative models for the free and forced motions 
of astronomical bodies (specially for the Earth) . Historically the perturbed Newtonian Euler equation for an elastic 
deformable Earth was apparently first derived by Jeffreys and Vicentefl]. Some different forms of it have appeared 
later in the literature [2]. Such a local treatment of global geodynamics has been pursued especially by Wahr[3], 
Schastok[4], and Dehant & Defraigne[5] to describe the nutation of the Earth. Clearly all of these investigations 
just mentioned arc fully within Newton's theory of gravity. The theory of elasticity is also of importance for the 
interpretation of data resulting from modern observational techniques such as Very-Long Baseline Interferometry 
(VLBI) , Lunar and Satellite Laser Ranging (LLR & LSR) , and all other kinds of observations where the positions of 
Earth-bound points should be described with high position. The normal modes (or quasi-normal modes) of the Earth 
or other astronomical bodies, such as white dwarfs or neutron stars is another field of important applications [6, 7, 8] 
as also the calculation of the time evolution of the (mass and current) multipole moments[9] of astronomical bodies. 

Extending the Newtonian theory of motion of elastic deformable bodies to include relativistic effects presents a 
new and improved basis for further discussions of the such problems. Recently Xu, Wu and Soffcl[10] developed 
such a general relativistic theory of elastic deformable astronomical bodies on the basis of the Damour-Soffcl-Xu 
(DSX) formalism as the foundation of modern general relativistic celestial mechanics at the first post-Newtonian 
approximation of Einstein's theory of gravity [11, 12, 13, 14]. In [10] we discussed the post-Newtonian perturbations 
of a uniformly rotating, stationary, and axisymmetric elastic body in a rotating Cartesian coordinate system. The 
general perturbations of such a configuration are treated within the Carter-Quintana formalism[15, 16]. A central 
result was the post-Newtonian dynamical equation for the displacement field in Cartesian coordinates representing 
the post-Newtonian version of the well-known Jeffreys- Vicente equation[2]. However, for practical applications, the 
common way of dealing with such perturbations is to go to spherical [17] instead of Cartesian coordinates and then 
to expand the relevant quantities in terms of scalar, vector and tensor spherical harmonics or so-called generalized 
spherical harmonics. 
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In this paper, we follow the route of our previous paper [10] to deduce the post-Newtonian perturbed local evo- 
lution equations and the perturbed Eulcrian equation for the displacement field of an elastic astronomical body in 
rotating spherical coordinates. To this end the Eulerian variation of Einstein's energy-momentum conservation law 
is performed. The Newtonian version of our results (neglecting all 1/c 2 terms) is in agreement with standard results 
from the literature (e.g., Ref.[17] after correction of a typographical mistake). Post-Newtonian junction conditions 
for the transition from one layer to another with different elastic properties that were not treated in [10] are also pre- 
sented here. Such general relativistic junction conditions have been discussed in a different context [18 , 19]. However, 
the application to a post-Newtonian displacement field presented here is new. Corresponding Newtonian junction 
conditions can e.g., be found in Wahr[3]. 

Symbols and notations are taken from the DSX papers [11, 12, 13, 14]: the space-time signature is taken as — + ++, 
spacetime indices go from to 3 and are denoted by Greek indices, while spatial indices (1 to 3) are denoted by Latin 
indices. We use Einstein's summation convention for both types of indices, whatever the position of repeated indices. 
We shall often abbreviate the order symbol 0(c~ n ) simply by 0(n). Local coordinates X a = (cT,X a ) are chosen in 
that reference system that moves with the body under consideration. The DSX scheme provides a description of the 
metric tensor G a p in a local system with two metric potentials (W, W a ). 

In Section II, the post-Newtonian perturbation equations of clastic astronomical bodies both in Cartesian [10] and 
spherical coordinates are presented. Special emphasis is given to the new derivation in spherical coordinates. In 
Section III, we discuss post-Newtonian junction and boundary conditions. In the last Section some conclusion can be 
found. 



II. POST-NEWTONIAN PERTURBATION EQUATIONS OF ELASTIC DEFORMABLE 

ASTRONOMICAL BODIES 



A. Perturbed equations in Rotating Cartesian coordinates 



The formalism starts by considering some isolated relaxed body that rotates uniformly with angular velocity 
about its symmetry axis with respect to some global non-rotating coordinate system. Deviations from such an 
equilibrium configuration and the action of tidal forces in a gravitational iV-body system are then described by 
means of perturbation theory. In that case one has to deal with at least three different coordinate systems: a global 
coordinate system x^ = (ct, x 1 ) like the Barycentric Celestial Reference System (BCRS) that extends to infinity and 
where the dynamics of the iV-body system can be formulated, a local 'non-rotating' system X a = (cT, X a ) like the 
Geocentric Celestial Reference System (GCRS) (e.g., kinematically non rotating or at least slowly rotating with an 
angular velocity of post-Newtonian order with respect to the global system) and finally some local coordinate system 
(X a = (cT,X a ) with T — T) whose spatial coordinates co- rotate uniformly with the equilibrium configuration. The 
post-Newtonian perturbation equations of an elastic deformable astronomical body in rotating Cartesian coordinates 
have been presented recently[10]. Here, some of the results are summarized. The components of the metric tensor in 
rotating Cartesian coordinates read 

— / 2W + V 2 \ 

Goo = -cxpf- y J +0(5), (2.1) 

— V a 4W„ 

G 0a - ^ + 0(5), (2.2) 

c c? 

( 2W\ 

G ab = ^exp — +0(4), (2.3) 



where 



and 



2WV 2 AW b V b V 4 

W = W+^-^- + ^ (2.4) 



W a = R ab (w b - \v b W^j . (2.5) 

A bar on top of some quantity indicates that it refers to the rotating coordinate system X a , otherwise it will refer 
to the 'non-rotating' local coordinate system X a . W and W a are the scalar and vector potential which describe the 
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metric in local 'non-rotating' coordinates. For the relaxed ground state of the body these potentials result entirely 
from the gravitational action of the body itself. For more details the reader is referred to Damour et al.[ll]. V b is 
the rotation velocity of the equilibrium configuration, R ab is a time-dependent rotation matrix (defined by (2.15)- 

(2.16) of [10]), V a = R ab V b and O = R bc VL c . fl c is the angular velocity with respect to ' non-rotating' coordinates 
(V a = e abc fl b X c , e abc is the completely antisymmetric Levi-Civita symbol of rank 3 with e 12 3 = +1). The perturbed 
energy balance equation after a first time integration reads (see (4.30) and (3.31) of [10]) 

5p = -V • (ps) - \ (pV°8 a + { V s a ). a + 2W. a s a p + ZpSW^j + 0(4) 

= - Pta s a -p*6 + 0(4), (2.6) 

where p is the energy density p* = p + p/c 2 the chemical potential per unit volume, p is the isotropic pressure and 
s a are the spatial components of the contravariant displacement field. The volume dilatation 6 is given by 

6 = s b M + 1 (V b s b + e bcd Tfs d V h + SW, c s c + 35W) + 0(4) . (2.7) 

Here, SW is the Eulerian variation of W in rotating coordinates. The perturbed Eulerian equation takes the form 
(Eq.(4.32) of [10]) 

= P* ( X + ^) (*° + ^abc^V) + P*®W G , a - P *S b a W G , b - P *(SW G ), a - p*S b W GM 

-{K®8 aP + 2psP a ).f, + 1 { P * [v a (V b M b )+W G , b s b V a - 2V b s b W G , a + (SW)^V a 
+8s b W [b , a] - 4{6W a ) >T ] + k (eW G , a - 9F a ) + (k(4H76 + vV/j) J + 0(4) , (2.8) 
where the post-Newtonian geopotcntial W G is given by 

_ — v 2 

W G = W+ — . 

The "dot" stands for the derivative with respect to the time variable T = T and 5 indicates the Eulerian variation. 
The elastic moduli k and p are the compression modulus and the shear modulus respectively; s ab is the shear-stress 
tensor (a complete representation of s ab is given in (4.26) of [10]). The term containing the shear-stress tensor reads 
(see also (4.33) of [10]) 

{2ps l3 a ) :[j = (2ps ba ), b + ^{-(4pWs ba + 2ps ca V c V b ), b + (2ps ab ). T V b 

+2p(2W, c s ac + e acb n c s bd V d ) } . (2.9) 

Equation (2.8) is the post-Newtonian Euler equation for the displacement field s, sometimes called the post-Newtonian 
Jeffreys- Vicente equation. 



B. Perturbation equations in rotating spherical coordinates 

1. Unperturbed and perturbed projection tensor in rotating spherical coordinates 

The importance to formulate the perturbation equations in rotating spherical coordinates was already stressed in the 
Introduction. Note that both the Eulerian and the Lagrangian variation of a tensor do not necessarily preserve the ten- 
sor character, whereas the difference between them, the Lie derivative = A — S), does. This implies that one cannot 
simply transform the variational equations from Cartesian to spherical coordinates. For that reason all derivations of 
the perturbed energy equation and the post-Newtonian Jeffreys- Vicente equation have to be repeated in spherical co- 
ordinates. These rotating spherical coordinates (r, 6, <f)) are defined by (X = r sin 9 cos (f>, Y = r sin 9 sin <f>, Z — r cos 9) . 

The metric tensor in such rotating spherical coordinates (X^ = (cT, r, 9, <fi) with T = T) takes the form 

( 2W + V 2 \ 

Goo - - cxp - \ + °( 6 ) - ( 2 - 10 ) 
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V c 4W C \ 

G 0a = G a a = D ac I =- + 0(5) , (2.11) 



<■ c 3 



Cr,,,, = D ul ,ox V ( H^j +0(4), (2.12) 



where U a = {dX a /dlt^ , W a = (dX a /dX b ^j w\ W = W, V 2 = V 2 and 

1 

D ab =\ r 2 | . (2.13) 
r 2 sin 2 6» 

The corresponding inverse matrix D ab satisfies D ab D bc = 5%. For the 3-dimension quantities V a and W a we define 

V a = D ab V h , W a = D ab W b . 

In the following we will choose ft = fle z , with £1 being constant, so that 

V r = V r = V 9 = Vg = 0, V^^n, % = n r 2 sin 2 6» (2.14) 

and V 2 = Q 2 r 2 sin 2 0. From these quantities the Christoffel symbols and the orthogonal projection tensor (projecting 
into the 3-space of an observer that moves with the corresponding material element) can be derived. In Carter and 
Quintana's formalism [15] a body is described by means of a bundle of time-like world- lines in a four-dimensional 
space-time manifold. A mapping into a three-dimensional manifold that is composed of the various material elements 
of the body identifies the various "particles" of the body. The projection of any tensor onto the local rest frame of 
matter is achieved with the orthogonal projection tensor 

7„ = G„ + ^. (2.15) 

This tensor acts as a positive metric tensor on the tangent subspace orthogonal to the flow vectors. The corresponding 
inverse metric tensor is given simply by its equivalent contravariant form, since it satisfies ^ v ^ vo =7^- is the 
four- velocity of some material element as tangent vector to its world line and j^uU" = 0. 

In General Relativity, the unperturbed and perturbed states of a body are considered as two configurations in 
separate four-dimensional space-time manifolds. Usually one starts with canonical coordinates in both manifolds 
and — > + Ax^ maps the coordinates of a material element in the reference state onto the coordinates of the 
same element in the perturbed state, where Ax^ is the position coordinate displacement in 4-dimension space-time. 
The quantities £ M = Ax^ are called the four-dimensional displacement field. The Lagrangian variation is the variation 
of the field in terms of a coordinate system which is itself dragged along by the displacement Ax M ; it is denoted by 
the symbol A. The Eulerian variation denoted by S is the variation taken at a fixed point in 4-dimcnsion space-time. 
The relation between these two kinds of variations is given by S = A — C^, where stands for the Lie derivative 
along the displacement field The displacement field, for obvious reasons, will be defined in rotating coordinates. 
It is taken as ^ = (0,£ a ) (see [10]). 

The Euler variation of G M „ (denoted by h^ v ) results from (2.10): 

( 2W + V 2 \ 

/loo - «?oo = — I 1 Z j + ( 6 ) , (2.16) 

h 0a = SG 0a = -D ac + 0(5) , (2.17) 

h ab = 5G ab = D ab (^^j + 0(4) . (2.18) 

Other important quantities are the Eulerian variations of the 4-velocity and the projection tensor. They take the 
form 

SU° = 1(SW + KC.t) + 0(3) , (2.19) 
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5U a = e,T [ 1 + ^ ) + 0(4) , (2.20) 



5f7 = -SW + 0(3), (2.21) 



c 



SU a - £ afc <^ e,T 1 + + - 2We,T ~ 4W 6 + V b (SW + V a C, T ) \ + 0(4) (2.22) 



and 



<57°o = -\v a C T + 0(4), (2.23) 

5 7 °6 = -D bc e.T + 0(3), (2.24) 
c 

<5 7 a o = ~-e,T + 0(3), (2.25) 
c 

^7% = ^C,TD bc V c + 0(4). (2.26) 

2. Lagrangian strain tensor and shear tensor 

For a perfect elastic body, the energy-momentum tensor of the relaxed state is of the form 

T af3 = pU a U +n a 0, (2.27) 

where p is the rest energy density and p the isotropic pressure in the reference state. The perturbed configuration 
changes the energy-momentum distribution and geometrical shape with time and might experience tidal forces from 
other astronomical objects. The Eulerian variation of the energy-momentum tensor can be expressed as 

5T af} = U a Uf}5p + p5(U a Up) + ■yapSp + pSj aP - 2ps a/3 , (2.28) 

where s a p is the shear tensor and p is the shear modulus. The symmetric trace-free shear tensor is defined as 

S[iu — &i±u 2®7/^ i (2.29) 

where the volume dilatation is given by O — e^. 
The Lagrangian strain tensor is defined by 

= \ = \ (l^ - 7^) , (2.30) 
where 7*^ is the unstrained value. Restricting ourselves to linear perturbations, the strain tensor is given by [15] 

= {Kp + 2£ (Q;/3) ) . (2.31) 

In the following we choose the displacement field in rotating spherical coordinates as £@ — (0, £ r , £ e , and we use 
the linear displacement 

D = D r e r + D e e g + D^e^ (2.32) 

with 

D r = C, D 6 = r£, , D 4 ' = rsm6^. (2.33) 

The explicit calculation of the Lagrangian strain tensor e^ u and the shear tensor s M „ is straightforward but cumber- 
some. Results for e^ u and to post-Newtonian accuracy are given in the Appendix. 



3. The post- Newtonian energy and Eulerian equations 



The Eulerian variation of the pressure Sp can be derived similarly as in our previous paper dealing with Cartesian 
coordinates [10]. In spherical coordinates it takes the form 



Sp = -p*D • VW G - kO + - J (iOW + fl 2 r sin 0(£>* + sin 9D r + cos 0D B )) + 0(4) . 
Our main results concern the perturbed local evolution equations 

For jU = one derives the perturbed energy equation in the form 

s P = -d • v P - P o - -Ue 



(2.34) 



(2.35) 



= - (D r p, T + -D 9 Pfi + -^—D*p^ - p*Q + 0(4) . 
\ r r sins' / 



(2.36) 



The perturbed Euler equations correspond to the case p = a. As an intermediate result by using the expression for 
V • D from the Appendix we get 



V 2 D = e r 

+ e e 

+ e<t> 



r l r z smv r z smO 



V 2 D e 



r 2 sin ( 



V 2 D - 



r 2 sin^ 9 r 2 
D4> 



r 2 sin 2 I 



2 cos9 



r 2 sin 9 

where e r , eg and are unit vectors in the r, 9 and (j) direction respectively 

1 



r 2 sin 2 ( 



(V-D), r = D^ r -1^ + V er 



1 



sin (9 ^ r r 2 sin<9 



+ -IT r - ^L» r + - cot 9D e r - cot 9D e 



(V-B).e = D: r9 + -D e ee + -^—D\ 
r ' rsinO 

+ -cot9D% K-D 6 , 

r ' r sin 

(V-D), = D^ + -D% + -^D 4 



r sin 2 9 ^ r ' e 



-D r , + - cot 6D e 

r ,v r 



r '"^ rsin9 

Tedious calculations then lead to the explicit equations for p = a = r, 9, <p respectively. 
The p = r equation reads 



= p* 



1 + ^fj (p r ; TT - 2ftsin6»L> T ) + p*QW G , r - p*{6W G ), 



y(D.vw G ), r -( K e), r -2(^ 



2WD r TT + -D e T W [e , r ] 



r sin 



+4(K6TF),r + 



-D* T W [4 ,. r] - 2Q.rsm9W G>r D* T ~ 4(<W~),t 
1 



+ nQW G ,r 



kVL 2 ( + ^<dr 2 sin 2 ' 



0(4) ; 



(2.37) 
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where 



2( A iS /3 r ) ;/3 = + Wr;/J 

= /u Q(V • D), r + (V 2 D) r ) + 2/v ^ r - lv • d) 
+^(zr, + rD? p -^) + 



PM 



r sin V f sin 6* ,r r ) 



r 2 r 2 sin 



5 1 f2 2 

+2Os r0jT - ^ sin 0(rD* T ) ir - -0 2 r sin (9 cos <9i/ r - — -D^ 

o ^ 

- -ft 2 sin 2 6>(L> r + D e cot 6») - 2 sin OD* 

3 ' 

-^, r fir sin 9 (p* T + D r n sin 9 + D VL cos 6j + n^Sl (D> ~ nD ^) } + °( 4 ) • 



(2.38) 



Two indices enclosed in parentheses imply symmetrization as in W( 0) &) = (W 0) (, + M / b ;a )/2, and two indices enclosed 
in square brackets imply anti-symmetrization as in W[ a ^ — (W a ,b — W^ a )/2- 

The \i = 9 equation takes the form 

= p*r ^1 + ^rj (V TTt - 2ficos0D$,) + p*eW G , e - p*(SW G ),e 
-p*(T> ■ VW G ),e - (k9), ( - 2(/w£) ;/3 + 1 {p* [2r#D* TT + 8^ T ^ e] 



+ - 



r sin 



:I>? r W M - 2Qrsin0W G , .D T - 4($W 9 ),t 



+ rceW G , e 



+A(kQW), 8 



-6r 2 sin 2 ( 
3 



■0(4), 



(2.39) 



where 



2(/xs / Y). /3 = 2/^A + V 



mQ(V- D), fl + r(V 2 D)^ + /v (IT, + r£ e r - £ fl ) 



\ r r 3 / rsin0 \ sin0 



+ Dt - D 4, cot 9 



r 2 r 2 sm 2 6l 6 



-i« 2 r£>^ + ^ 2 r 2 sin0cos0£> r r - ift 2 rsin 2 0£> r fl - %tt 2 rcos9D* 



-2Q 2 r sin cos 9D r + -fiV sin 2 6>L> 9 - 2ft 2 r cos 2 9D e 
- ^n i6 flr sin 9 (p% + fl(D r sin 9 + D 6 cos 0)) 



(2.40) 



Finally, the p — 4> equation reads explicitly 

2WV 



= p*rsin0 1 + 



(d tt + 2ft sin 6(D r T + cot 9D e T )^ 
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where 



+P*QW Gt<t> - p*(6W G ) t<j> - P *(D ■ VW G )^ - (k&)^ - 2(^) ij8 

{/>* \2rsin6W G D% T + VLr 2 sin 2 6>(D, T ■ VVF G ) - 2VLr sm6W G:4> D% 

+%D r tT W [rM + ^D e T W [eM + nr 2 sin 2 0{5W), T - 4(JW^),t 



+k(6VF g , - 6, T 0r 2 sin 2 0) + 4(k9TT), + 



kO 2 I 



-0r 2 sin 2 i 
3 



+ 0(4), 



2(mA) ;/3 = + 2/i«^ ;/3 



= AiQ(V-D) :0 +rsin6»(V 2 D) ^ + /i >r .rsin0 

^ . ^ D*cote\ 

+/i^sm9 H ^— 

\ r rsmfr 1 r / 



rsinf ' r / 



+2fi, 
1 

+ - 
1 



rsinO r r 3 



2/i 



2s^W. r + 2 -^W. e + 4%-W 4 
r 2 r 2 sm 2 d 



: flr 2 sin 2 6»(V • D), T 



+ §fiV sin 3 0(V 2 D) + -fir sin 0£>^ t 
+20r sin 0(£T T sin 6» + D e T cos 6>) - ^Q 2 r sin 6 (D^ sin (9 + D ^ cos ( 
+ io 2 rsin6» (rsin 2 0D+ - D 4, + 2cos6D^ g - D* J 



+/i !r Or 2 sin 2 
+[ififlr sin 2 



Or sin 



r sin 6 



Dt- — )+ IFr - flD\ 



\ r r sm r 



+ |/i,*fir sin 6> (V T + 0(L> r sin 8 + D 9 cos (?)) | + 0(4) . 



(2.41) 



(2.42) 



Eqs.(2.37), (2.39), (2.41) together with (2.36) are the desired dynamical equations for the displacement field. They 
are valid up to terms of order 1 /c 4 and second order in the displacements field itself. The Newtonian limit of our 
results (neglecting all 1/c 2 terms) agrees with standard results from textbooks (e.g., the ones from [17] after correction 
of a typographical mistake). 



III. POST-NEWTONIAN JUNCTION CONDITIONS 



In most cases such post-Newtonian dynamical equations of elastic dcformable bodies will be applied to astronomical 
bodies composed of different layers. If the body has several layers as e.g., the Earth that shows a solid inner core, a 
fluid outer core, a mantle and a thin crust, we have to consider corresponding junction conditions at the interface of 
two adjacent layers. First we consider such junction conditions in Cartesian coordinates. For practical applications 
they are then also formulated in spherical coordinates where they can be compared with well-known Newtonian 
results [3]. Junctions conditions arc formulated in rotating coordinates and to have a well defined stationary and 
axisymmctric ground state it is assumed that all layers rotate with the same angular velocity f2, i.e. there is no 
relative motion between two layers in the ground state. The behavior of individual physical quantities and their 
corresponding perturbed quantities on the interface will be studied at first. 



9 



The gravitational potentials W and W a are physical quantities in the non-rotating ground state. For a isolated 
body, W and W a can be obtained as a solution of Eqs.(2.4), (2.5) in [10]. They are inhomogeneous D'Alembert's and 
Poisson differential equation respectively. Although the sources £ and T a may be discontinuous across any interface, 
W and W a (solutions of the equations) are continue on the interface. Since we do not consider the shock wave, the 
surface mass-density and surface current mass-density do not exist anywhere. Therefore derivatives of W and W a have 
to be continue on the interface also. W and W in rotating coordinates differ from (W, W a ) simply by a nonsingular 
coordinate transformation (Eq.(2.4) and (2.5)). So W, W and their derivative continuous on any interface as well. 

To discuss the continuity of SW and SW we have to consider SW and SW a in non-rotating coordinates before 
hand, since only in non- rotating coordinates SW and 5W a can be obtained from perturbed field equations (Eq.(4.17) 
and (4.18) of Ref.[10]) 



i d 2 \ 

c 2 dT 2 J 



SW = -4nGST + 0(4) (3.1) 



V 2 SW a = -AnGST a + 0(2) (3.2) 



where 

ST = -p* a s a -p*e 



-2, 



2p*(5W + 2V s a ) - 2{p* a s a + p*0){W + V ) 

L. L 

2p*W Gtb s b - 3k9] + 0(4) , (3.3) 



5T a = R 



ba 



p*s b -V"(p* c s c + p*&) +0(2). (3.4) 



SW and 5W a are solutions of Eqs.(3.1), (3.2). They are continuous as the same discussion on W and W a , but 
their derivatives arc in different cases. The Eulerian variation of the surface mass-density and Eulerian variation of 
the surface current mass-density do exist on the interface because both ST, and ST a are dependent on the spatial 
coordinate derivatives of p (to see Eq.(3.3), (3.4)). Then on interface ST, and ST a are divergence. Therefore (5W) ya 
and {8W a )^b are finite on the interface but not necessary to be continuous across the interface. Through the coordinate 
transformation and neglecting all of higher order terms, we get the representation of SW and SW by means of SW 
and SW a (see (4.22) and (4.23) of [10]) 

SW = SW+^ (2SWV 2 - 4SW b V b ) + 0(4) , (3.5) 
8W a = R ab [5W b - \v b 8W ) + 0(2) . (3.6) 



2 



The inverse transformation takes the form: 



SW = SW + ^V a SW a + 0(4), (3.7) 

SW a = R ba (8W h + ^V b 5W^j + 0(2) . (3.8) 

The transformation formulae is nonsingular, then the behavior of SW and SW a are similar to SW and 5W a , i.e. 
SW and SW are continuous across any interface. SW. a and (SW ),b are finite on the interface, but not necessary 
continuous. The first and second time derivatives of SW and SW are continuous everywhere including in the interface 
since the shock wave does not exist at anytime. Later when we consider junction condition, we can drop the continuous 
terms on both sides of the interface. 

For the displacement field s we shall take a similar physical consideration as in [3] . The field s is continuous across 
any solid-solid interface. Its normal component n • s is continuous across any interface. But its tangent component 
maybe not continuous for solid-liquid interface, since the tangent interaction between solid and liquid is close to zero 
in the absence of viscous forces. Also s a and s° are finite across any interface, p, k, s a b and <d are finite, but on the 
different side of a interface they maybe different (not continuous). Therefore p^ a , Sbc,a and iCl are not necessary 
finite on the interface. As we mentioned before all of layers rotate with the same angular velocity f2, so that V a and 

V a (V a = R ab V b ) are continuous across the interface. V h a , V a and V 2 = V a are continuous as well. 

In Eq.(3.1) the DAlembertian can be related to the non-rotating coordinates. A transformation to rotating coor- 
dinates yields 

i d 2 -2 2V a d 2 v a v b d 2 i d 2 

~ V + „2 ^ a ^7TT „2 ^h^a „2 „=2 ' ^' 9 ) 



c 2 dT 2 c 2 dx a df c 2 d x h dx a c2 dr 2 
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where we have used the relation dV a /dX a = 0. Substituting Eq.(3.9) and (3.7) into Eq.(3.1), we get 



d 



d 



■SW 



-zV SW^ 

c 2 1 c z dX 



dx 

= -AirGSH + 0(4) , 
We also can rewrite SY, (to see Eqs.(4.27), (4.31) and (4.40) of [10]) 



V V dSW 
c 2 dt 



T sw 



,TT 



(3.10) 



SY 



-V • (pa) 



1 



3 P V a s a - (ps a ) ia - 3 P W ib s b - P 5W 



-2V • (ps{W + V 2 )) + \pV 2 .bS b - 3ns\ b 



(3.11) 



From now on for convenience we omitted "bar" in V as we did in [10]. Substituting Eq.(3.11) into Eq.(3.10), we get 



V • A + B = 0(4) 



(3.12) 



with 



VSW - AnGp*s 
1 



+ - 



2VSW, Y + 4V(F SW ) - V(V • VSW) - %kG P *s(W + V ) 



and 



B = 



(l27rGps • 



tt2 



V - UvGps ■ VW - 47rGp5W + 2nGps ■ \7V - UirGne - SW ™ 



As the discussion on the boundary condition problem in classical physics, we now integrate Eq.(3.12) over an infmites- 
imally small volume AV such that an interface intersects this volume. We choose AV as a cylinder. The surface (S) 
of AV enclose the interface of two different layers with the shape of a circular drum of a radius r and a depth h, 
where h <C r and r is so small that the portion of the interface contained in (S) can be taken as flat. When h — > 
the integration divide into two parts: the first part of Eq.(3.12) becomes a surface integration by means of Gauss' 
theorem, the second part tends to zero with vanishing AV because the terms in B are all continuous or finite. Hence, 



n • A|i ayor i — n • A|i ayer 2 . 



(3.13) 



i.e. 



' ' b 5W b ) - V(V • VSW) - 8nGp*s(W + V 2 ) }+ 0(4) 



n- jv<SW-47rGp*s+ ^ 4V(y 
continuous across any interface , 
here we have dropped the term 2V(SW),t since it is continuous across interface. 



(3.14) 



Substituting Eq.(3.8) and (3.4) into (3.2) and considering V = V 2 and the rotation matrix R ba dependent on time 
only, we get 



d 



®( 5 W a + l -V a SW) ^Gs b pV a 

dx V 2 



inG Ls a + ps"^) =0(2). 



Performing a similar integration over an infinitesimal volume as before, we find that 



n • 



V(SW a + ^V a 5W) - 4ttGp*sV° 



0(2) continuous across any interface , 



(3.15) 



(3.16) 



where p has been substituted by p* since this equation is valid only to 0(2). Eq.(3.16) can be written by means of 
the form of parallel vector, i.e. 



n ■ 



V(<5W + \SW) - 4ttG P *sV 



+ 0(2) continuous across any interface . 



(3.17) 
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Finally a similar integration of Eq.(2.8) over an infinitesimal cylindrical volume on the interface leads to 

{Aba,b + B a )dV = , 



Ay 



where 



1 /. ^ 

72 



!,„, = - iA-)6 ab - 2ps ba + ^ ( 4pWs ba + 2fis ca V C V b + 4 K W@S ab + kV"V 's d c 6 ab 



) 



(3.18) 



(3.19) 



B a = p* 1 + 



2W ( 



(s a 



+ 2e abc tfs c ) + p*QW G , a - p*s h a W G M - p*(SW G ), a - p*s b W GM 



+ 4, { P * [v a (V b s b ) + W G Ms b V a - 2V b s b W G , a + (5W)^V a + 8s b W [b , a] - 4(SW a ); 
-(2ps ab )^V b - 2p(2W, c s ac + e acb n c s bd V d ) + K{®W G , a - 67°)} . 
All of terms in B a are finite or continuous as we mentioned before. In terms of median method, we have 



B a dV = B a irr 2 h. 



AV 



The first term of Eq.(3.18) can be deduced as an surface integration 

/ A ba , b dV = irr 2 {n\A ab - n\A ab ) . 

J AV 

Cancelling 7rr 2 and neglecting the higher-order term B a h, we get 



(3.20) 



(3.21) 



(3.22) 



+ 2ps ab - 4r 



4pWs ba + 2ps ca V c V b + A K WQ8 ab + kVVV, c <U } + 0(4) 



continuous across any interface . 



(3.23) 



We should point out that the Newtonian part of Eq.(3.23) n05 ab + 2ps ab is just the Newtonian stress tensor T ab in 
[3], therefore Eq.(3.23) is an extended PN version of Wahr's. 



The PN junction conditions in Cartesian coordinates are summarized as follows: 

s continuous across any solid — solid interface 
s • n continuous across any interface 
SW continuous across any interface 
5W a continuous across any interface 
1 
6 

continuous across any interface 



n - jvcM'. - a.tC'p 's 



4V (VW fc ) - V(V • V(5W)) - 8ttG P *s(W + V 2 ) j + 0(4) 



1, 



V((5W + -V5W) - 4ttGp*sV 



0(2) continuous across any interface 



j«6<5 ai 



b + 2ps ab - 



ApWs ab + 2ps ac V V + 4 K WeS ab + K S ab V V s d , c } + 0(4) 



continuous across any interface . 



(3.24) 
(3.25) 
(3.26) 
(3.27) 



(3.28) 
(3.29) 

(3.30) 



For the outer surface boundary conditions wc only need simply take p = n = p = s = outside the elastic body, i.e. 

(3.31) 



«5VF| in = 6W\ out 

sw a \ in = sw a \ out 



(3.32) 



n-|v5W-47rGp*s+^ 4V (vW 6 ) - V(V • V(SW)) - 8irGp*s(W + V" 2 ) J 
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4V (VW fc ) - V(V • V(<W))] I lout + 0(4) 
n • | V (W + 1 VtfW^ - 4nGp*sV^ \ m = n • V + ^ VW) 



0(2) 



+ 2^s ah - 



4 M iy Sa6 + 2 f u ac V c V b + A K WeS ab + K S ab V d V c s d , 



0(4). 



(3.33) 
(3.34) 

(3.35) 



For most applications the use of spherical coordinates will be advantageous. Corresponding boundary and junction 
conditions cab be derived similarly to the case of Cartesian coordinates. They read 



D continuous across any solid — solid interface 
D • n continuous across any interface 
SW continuous across any interface 

SW a continuous across any interface (a = r,6,(j> respectively) 

n • | V6W - 4wGp*T) + 4V(ftr 2 sin 2 0SW*) - 8nGp*T>(W + fiV sin 2 9) 

— e<j,Cl 2 r sin 9 (8W) t< /, j + 0(4) continuous across any interface 

0(2) continuous across any interface 



n 



V (SW + ^VSWj - 4ttGp*DV 



n b I K QD ab + 2ps ab - -2 \4nWs ab + 2ps atj> ti 2 D H + 4 K W@D ab + nD ab (fi 2 r sinflD*, 



+ -D- V(Q 2 r 2 sin 2 9) 



continuous across any interface . 



(3.36) 
(3.37) 
(3.38) 
(3.39) 



(3.40) 
(3.41) 

(3.42) 



s ab in Eq.(3.42) is the shear tensor in spherical coordinates, which is shown in the Appendix Eqs.(A.10)-(A.15) 

For the outer surface boundary we take /i = k = /j = D = outside the elastic body as before. When we neglect all 
of 1/c 2 terms and let V = O(non- rotating), our formulae (Eq.(3.40) and (3.42)) agree with Wahr's results (to see Sec. 
Ill of Ref.[3] ), i.e. our work is an extension of the Newtonian version to rotating PN version. As for our formulae 
(Eq.(3.41)), since it is a purely PN junction condition, there is no Newtonian corresponding to be compared. 



IV. CONCLUSION 



In this paper we present the perturbation equations for the dynamical behavior of astronomical elastic bodies in the 
first post-Newtonian approximation of Einstein theory of gravity in rotating spherical coordinates. In comparison with 
our previous in Cartesian coordinates [10], these equations in spherical coordinates are more useful for applications 
since usually all relevant equations are expended in terms of scalar, vector and tensor spherical harmonics. The 
equations and relations such as junction conditions can e.g., be applied to problems of geodynamics (e.g., for the 
theory of nutation) or seismology of compact stars (e.g. for the problem of the normal modes of astronomical bodies). 
Also the formulation of post-Newtonian junction conditions at the interface of adjacent layers of different elastic 
properties is presented here for the first time. A comparison reveals that the Newtonian limit agrees with well-known 
results from the literature. Here it should be emphasized that the junction conditions can be written very generally 
with ordinary Euclidean 3-vectors and 3-tensors so that they can be formulated for a broad class of coordinate systems. 

Our perturbation equations together with junctions conditions can in principle be solved if the internal quantities 
of state (density, pressure etc.) and elastic moduli are given e.g., by some Earth's model[20]. However, these 
equations describing the free and forced motions of the body are complicated partial differential equations. For that 
reason usually an expansion of relevant quantities in terms of spherical harmonics turns these equations into a set of 
coupled ordinary differential equations. A different natural basis for such an expansion is provided by the so-called 
generalized spherical harmonics [21, 22] that was employed by Wahr and other authors. By using the generalized 
spherical harmonics to expand all functions (displacement vector, incremental Eulerian gravitational potential energy, 
incremental elastic stress tensor, applied force (tidal force et al.)), the original partial differential equations and 
boundary conditions are transformed into a set of ordinary differential equations and scalar boundary conditions 
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for the unknowns functions. For a spherical and non-rotating ground state these ordinary differential equations are 
uncoupled, for an oblate, rotating body they are coupled, the coupling parameter being given by the dimensionlcss 
oblatencss of the body. Such expansions together with applications for geodynamics will be published separately. 
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APPENDIX: EXPLICIT RESULTS FOR THE STRAIN AND SHEAR TENSOR 



From (2.31) one derives 

eoo - 0(4) : 



eor = em = eo<f> = 0(5) , 

e rr = D r r + ^ (SW + D • VW + 2WD^ + 0(4) : 



1 



1 



1 



2W 



2 • 2 

= r sin i 



■ sin 9 



1 

r 



D r + -D e cote + \ 



2W 



-D% + -D r + —\6W+ D r + D • VW + D 



0(4), 



SW + D • VW 



+ — ?-r (2W + V 2 )D' t ' lh + nrsm9D 4 ' T + -(W + V 2 )(D r + D 6 cot 6) 
rsmO ' v ' r 



± iD:9 + rD° r -D°) ( 1 + 2 J)\ +0(4), 



0(4), 



(A.l) 
(A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 



1 / D 



= rsm( 



4> 



2 V r sin t 



-D^ + D+ 



W 



r sin( 



n4> i i 

+(£>£ - —){W + -V 2 ) + —Or sin 6D r T 



+ 0(4), 



1 / D e 



= r smt 



2 \ rsin( 



+ -D\ 

r r ' 



W 



r sin( 



+ 0(4) 



(A.7) 



(A.; 



where 



D • VW = W r D r + -W D e + -^—Wa.D 4 ' . 

r ' rsmV ' 



The volume dilatation 0, which is related with the expansion rate 9 by 9 = Q^U^, reads 



9 = e% = G^e^ = G ab e ab + 0(4) 



V-D + 4j [35W + 3B- VW + Slrsin9D' t ; T + )-V 2 t c ) +0(4), 



(A.9) 



where 



V D = D r r + -D% + 



-— D*+-D r + - cot 6D» , 
r sin ft ' v r r 



1. 



VX = Q 2 r sin 2 flf + fiV sin 9 cos 6? = sin 2 8(D r + cot 9D a ) . 
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From this it is easy to derive the symmetric trace- free shear tensor s^. One finds 



and 



*oo = 0(4), 

SOr = S 0e = S O = 0(5) , 
Srr = D^-^V-D + l 



2W D 



-- (fir sm 8 D% + ti 2 r sin 2 6 (D r + cot 6 D & 
3 



0(4). 



see = r 



-D r 



1. 



-V D 



2W I -D g + -D r 

r ' r 



:6 



i ^rsm6»L> T + fi 2 rsm 2 6i(D r + cot6lL> e )) | + 0(4) 



2 • 2 , 

r sin ( 



1 



1 



r sin# 
(2W + V 2 ) 



■ / ) 1 + il) r + -D 9 cot 6 — -V • D 

r r 3 



^-r^ + V + V cot 0-^9 



+fi 2 r sin 2 (9 (D r + cot 8D e ) + -fir sin 6»Z^ T 



--r2 2 rsin 2 0(£> r + cot 0D e ) 



+ 0(4) , 



(A.10) 
(A.11) 



(A.12) 



(A.13) 



(A.14) 



s r e = e r e , 



5r0 



se<p = ee0 . 



(A.15) 



These are the post-Newtonian components of the Lagrangian strain and shear tensor in spherical coordinates. If all 
post-Newtonian terms are neglected (i.e., for c — > oo), they reduce to the well known Newtonian results that can be 
found in the standard literature. 
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